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Abstract 

A broad range of single field models of inflation are analyzed in light of all relevant 
recent cosmological data, checking whether they can lead to the formation of long-lived 
Primordial Black Holes (PBHs). To that end we calculate the spectral index of the power 
spectrum of primordial perturbations as well as its first and second derivatives. PBH 
formation is possible only if the spectral index increases significantly at small scales, i.e. 
large wave number k. Since current data indicate that the first derivative 05- of the 
spectral index ns{ko) is negative at the pivot scale /cq, PBH formation is only possible 
in the presence of a sizable and positive second derivative ("running of the running") 
/3s- Among the three small-field and five large-field models we analyze, only one small- 
field model, the "running mass" model, allows PBH formation, for a narrow range of 
parameters. We also note that none of the models we analyze can accord for a large and 
negative value of as, which is weakly preferred by current data. 
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1 Introduction 



An epoch of accelerated expansion in the early universe, inflation, dynamically resolves many 
cosmological puzzles of the hot big bang model, such as homogeneity, isotropy and flatness of 
the Universe pQ. On the other hand, the generation of a spectrum of primordial fluctuations 
in the early universe is a crucial ingredient of all inflationary models. These fluctuations can 
explain the generation of all (classical) inhomogeneities that can be seen in our universe, from 
the Cosmic Microwave Background (CMB) anisotropies to the Large Scale Structure (LSS) in 
the form of (clusters of) galaxies [2]. The inflationary paradigm therefore reconciles Big Bang 
cosmology with the appearance of an inhomogeneous universe. In addition, many models of 
inflation |3j make accurate predictions allowing for observations of ever increasing variety and 
quality to discriminate between the various candidates. 

One such prediction concerns the possible formation of Primordial Black Holes (PBHs). 
Indeed, it was realized already some time ago that the spectrum of primordial fluctuations 
could lead to the production of PBHs [H [5]. For this generation mechanism to be efficient, one 
typically needs a "blue" spectrum [HI E] . This means that the amplitude of density fluctuations 
must be much higher at the small length scales relevant for PBH formation than at the much 
larger length scales probed by data on the Cosmic Microwave Background (CMB) and Large 
Scale Structure (LSS). This discrepancy of scales makes PBH formation also a unique probe of 
cosmological inflation; in particular, the constraint that not too many PBHs should have been 
produced has been used to limit the power spectrum at small length scales, which in turn allows 
to constrain models of inflation [S^iOj. For single-fleld inflation models, the relevant parameter 
space is deflned by the scalar spectral index ns; the ratio of tensor to scalar fluctuations r; the 
running of the scalar spectral index as', and we will see that the "running of running of the 
spectral index", is also important for PBH formation. 

The goal of this paper is to make use of the recent observational bounds on these parameters 
derived from the combined CMB data from the Wilkinson Microwave Anisotropy Probe seven 
year data (WMAP7) [lO] and the South Pole Telescope (SPT) [TT], as well as data on the Baryon 
Acoustic Oscillations (BAO) [12], the Hubble constant Hq [13], and clusters of galaxies [H]. 
We will investigate a wide range of inflationary models, checking whether they still can give rise 
to signiflcant PBH formation given these constraints. In so doing, we also check whether these 
models can account for a sizably negative running of the spectral index, as (weakly) favored 
by current data. We focus on models where the cosmic expansion was driven by a single, self- 
interacting scalar inflaton fleld 0. Moreover, we only consider models with simple potentials, 
which have been suggested for reasons not related to PBH formation. 

This paper is organized as follows: In section 2 we present a brief review of the Press- 
Schechter formalism describing PBH formation. In section 3 we flrst summarize the current 
bounds on the observational parameters, and review their calculation from the potential of the 
inflaton fleld. We then systematically analyze three small-fleld and flve large-fleld models. In 
section 4 we present our conclusions. 
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2 Primordial Black Holes Formation 



PBHs are black holes that result from the collapse of density fluctuations flGl [T7j. They are 
very sensitive cosmological probes for physics phenomena occurring in the early universe. They 
could be formed by many different mechanisms, e.g., from initial density inhomogeneities [TB] , 
a softening of the equation of state |TB] {e.g. in phase transition or during the preheating period 
after inflation), collapse of cosmic string loops [I9], bubble collisions [20], collapse of domain 
walls [21], etc. 

The idea that large amplitude matter overdensities in the early universe could have collapsed 
through self-gravity to form PBHs was first studied by Zel'dovich and Novikov [16], and then by 
Hawking ^/Tj. This theory suggests that large amplitude inhomogeneities in the early universe 
overcome internal pressure forces and collapse to form black holes. A lower threshold for the 
amplitude of such homogeneities 5th = (5p/p)th, was first provided by Carr [5j, giving 5th ~ 1/3 
at the time of radiation domination (RD){^The probability of PBH formation is a useful tool 
to constrain the mean amplitude of inhomogeneities on scales which cannot be probed by any 
other method. Since PBHs behave like matter, their contribution to the energy density increases 
with time during the RD epoch. For this reason, the PBHs formed considerably before the end 
of RD are the most relevant to cosmology. We have focused in our study [7] on these kind of 
PBHs and we also considered the standard case of PBHs formation, which applies to scales 
which have left the horizon at the end of inflationj^ We only consider Gaussian and spherically 
symmetric perturbations and we assume that the mass of the PBH formed is proportional to 
the mass of the horizon mass at horizon entry, Mpbh = T^pnj^ 

Generally, for Gaussian primordial fluctuations, the probability density P{S] R), where 6 is 
the density contrast averaged over a sphere of radius R, is given by 

Here, the dispersion (mass variance) (Ts{R) is computed using a Gaussian window function 
W{kR) = exp ' 



al{R) = W\kR)Vs{k)^, (2.2) 

where Vs{k) is the power spectrum of 6 which is related to the power spectrum of curvature 
perturbations on comoving hypersurfaces as follows [2]: 

^.(M) = #^(^Vp..W. (2,3) 



(5 + 3w)2 \aH_ 

Therefore the probability /(> M) that a PBH with mass > M is formed on a scale R when 



^Niemeyer and Jedamzik |22) carried out numerical simulations and found the threshold for PBH formation 
to be 0.7. We have shown [7] that PBHs abundance is sensitive to the value of i5th- 

^It has been shown [5] that PBHs can also form on scales which never leave the horizon during inflation, and 
therefore never become classical. We do not consider this contribution. 

•^A simple analytical calculation suggest that 7 ~ w^^^ ~ 0.2 during the radiation era [Sj. 
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that scale reenters the Hubble radius, is given bj|^ 

f{>M) = 2-f P{6;M{R))d5 . (2.4) 

In order to complete the calculation we just need to relate the PBH mass M to the comoving 
smoothing scale R when the scale enters the horizon, R = {aH)~^. It is straightforward to show 
that in the RD era 

« 5.54xlO-Vi(^V'Vj?^V", (2.5) 



1 Mpc ' V 1 g / ^3.36. 

where is the effective number of relativistic degrees of freedom which is expected to be of 
order 100 in the early universe. 

The power spectrum of primordial scalar fluctuations at scales A; ~ A;^ is typically parame- 
terized as a power-law with power ns, V-jidk) = 'PndkRj^k/ kR)"^^^^^'^ , with = 1/R. With 
this ansatz, the variance of the primordial density field at horizon crossing is given by 

^KR) = ^J^, VnAkR)n{ns{R) + 3)/2] , (2.6) 

for ns{R) > -3. 

In order to relate the scales relevant for PBH formation to the scales probed by data on the 
CMB and LSS, we parameterize the power spectrum as 

VTzSkR) = VnXk,){kR/kor^^^-^ . (2.7) 

It is important to distinguish between ns{R) and n{R) at this point. ns{R) describes the slope 
of the power spectrum at scales k k^ = 1/R, whereas n{R) fixes the normalization of the 
spectrum at kji ^ feg- The two powers are identical if the spectral index is strictly constant, 
i.e. if neither ns nor n depend on R. However, in this case CMB data imply pLOj that n = ns 
is close to unity. Eqs.(2.6) and (2.7) then give a very small variance, leading to essentially no 
PBH formation. 

Significant PBH formation can therefore only occur in scenarios with running spectral index 
[23] . We parameterize the scale dependence of n as [2 



n{R) = nsiko)-^as\n{koR) + ^Ps\n''{koR) + ...; (2.8) 

recall that we are interested in i? <^ 1/ko, i.e. \n{koR) < 0. The parameters as and f3s denote 
the running of the effective spectral index ns and the running of the running, respectively: 



nsiko) 
as{ko) 

Psiko) 



dink 

dns 



k=ko 



dink 

d^ns 



dln'k 



k=ko 

(2.9) 

k=ko 



*We will show that in practice P{S; R) is such a rapidly decreasing function of S above 6th that the upper 
cutoff Scut is not important. 
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Eq.(2.8) illustrates the difference between n{R) and ns{R)- The latter has an expansion similar 



to eq.(2.8), but with the usual Taylor-expansion coefficients, 1 in front of as and 1/2 in front 
of /3s. One therefore has 



ns{R) 



n{R) 



^ as In (fcoi?) + ^ /3s In^ {k^R) + • • • • 



(2.10) 



Setting ns{ko) = 1 for simphcity, eq.(2.10) implies ns{R) = 2n{R) — 1 for j3s = 0, and 
ns{R) = 3n{R) — 2 for as = 0. We will compute the variance cr{R), and hence the PBH 
fraction, for the former relation. 
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Figure 1: Fraction of the energy density of the universe collapsing into PBHs as a function of 
the PBH mass. 

The result of this calculation is shown in figure [TJ Here we have fixed 7 = 0.2, and show 
results for the threshold ^th = 1/3- Due to Hawking radiation PBHs contributing to Dark 
Matter (DM) today must have Mpbh ^ 10^^ g; at this mass, they saturate the DM rehc density 
if / ~ 5 X 10~^^ [26]. Figure [1] shows that this requires n{R) ^ 1.37 for 5th = 0.3. In order to get 
long-lived PBHs the amplitude of the perturbations at PBH scales must therefore exceed that 
at CMB scales by a factor 10^-10^. Current data favor a negative or at best slightly positive 
value of ^5 at the CMB pivot scale, as well as a spectral index at the pivot scale somewhat 
below 1. The first two terms in eq.(2.8) can thus not lead to PBH formation; so in the next 
section we study single-field inflation models by also considering the "running of the running 
of the spectral index" for PBHs formation. When possible, we will compare this to an exact 
calculation of the power of density perturbation at scales relevant for PBH formation. 
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3 Inflation Models and Primordial Black Holes 



Most models of inflation predict an approximately scale-free spectrum with a spectral index 
ns (as well as n) close to the scale-invariant (Harrison-Zel'dovich) case ns = n = \. As shown, 
a significant number of long-lived PBHs can only be produced for n > 1 (a "blue spectrum"), 
since these values lead to more power on small scales |i6j. Observational limits (both from 
Hawking radiation and the fact that PBHs must not overdose the universe) strongly constrain 
ns [8]. This, therefore, yields a constraint on inflationary models that is independent of the 
cosmological constraints from the CMB and LSS. However, we will see below that in most 
simple models of inflation constraints on model parameters derived from the latter are far more 
stringent than the PBH constraint, to the point of making the formation of long-lived PBHs 
impossible. 

A given inflation model can be described to lowest order in slow roll by three independent 
parameters: the normalization of the curvature perturbation spectrum some scale, the 

tensor-to-scalar ratio r, and the spectral index ns- If we wish to include higher-order effects, 
we have the forth and fifth parameters describing the running as-, and the running of running 
of the spectral index Ps- 

Observational bounds on V-ji^., r, ns and as at the pivot k^imt = 0.015 Mpc~^, where ns 
and as are essentially uncorrelated, are reported in |Tl] as follows^ 



ns(A;pivot) = 0.9751 ±0.0110, 

as{kp„ot) = -0.017 ±0.012, 

VT^^{ko) = (2.33 ± 0.092) X 10"^ 

r < 0.17(95%CL). (3.1) 



By requiring ns{k) G [0.9531, 0.9971] and as G [-0.041, 0.007] for all k G [10-^ 10] Mpc" 



down to the Lyman-a range), we find from eq.(2.8) that values of (3s up to 0.017 are allowed. 



The ranges of and a s are simply in 2 cr range of ( |3.1 ) , and the range of k encompasses 



all cosmologically relevant scales. Here we used the error bars derived from the analysis of 
the "SPT±WMAP7±BAO±i/o± C^^isters" data setQ This upper bound on Ps is lower than 
the bound found from the "WMAP+i^o+BAO" data set [7j, since the older data set allowed 
somewhat larger values of as- Note that the estimates for spectral index and its running are 
highly correlated for the typical COBE scale, fco = 0.002 Mpc~ 



Eq.(2.8) shows that for = we only need Psik^) ~ 0.0015 in order to generate sufficiently 
large density perturbations to allow formation of 10^^ g PBHs. Even if we set as{ko) equal to 
its central value, asiko) = —0.017 , we only need f3s{ko) ~ 0.0028. In this model-independent 
analysis including the running of the running of the spectral index thus easily allows to ac- 



^The amplitude of the primordial scalar fluctuations is reported at the "COBE" scale ko = 0.002 Mpc^^. 

^We ignore possible tensor modes, which is appropriate for small-field models. Allowing a sizable contribution 
from tensor modes changes the mean value of as [101; but unfortunately the pivot scale where the spectral 
index and its running are uncorrelated is not reported in running+tensor model. In the SPT data, the inflation 
parameters are not reported in running+tensor model. So although in the large-field models, tensor modes are 
not negligible we will assume that the upper bound of f3s is the same as in small-field models. Note that the 
precise value of the upper bound on derived here is not important for our analysis, since it is in any case 
well above the lower bound needed for successful PBH formation. 
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commodate PBH formation in scenarios that reproduce all current cosmological observations 
at large scales. 

Of course, this kind of model-independent analysis does not show whether simple, rea- 
sonably well-motivated inflationary models exist that can generate a sufficiently large Ps- In 
the following we study different models of inflation and check whether they can lead to PBHs 
formation. As a by-product, we also check whether these models can accommodate a sizably 
negative value of as, as indicated by current data. 



In order to calculate the spectral parameters ns, as and (3s deflned in eqs.(2.9), we need 
the flrst four slow-roll parameters, deflned as [2]: 





Ml fV'\ 


e 


2 \V) 


7] 


V" 


e 









A V"" 

2M^e— . (3.2) 

Here V{(j)) is the inflaton potential, Mp is the reduced Planck mass, and primes denote deriva- 
tives with respect to (p. Note that e is positive by deflnition whereas in spite of the square, 
can be either positive or negative. The square in and cube in are to indicate that they 
are second and third-order in the slow-roll expansion, respectively. These parameters must be 
less than one for the slow-roll expansion to be valid. All these parameters are in general scale- 
dependent, i.e. they have to be evaluated at the value of that the inflaton fleld had when 
the scale k crossed out of the horizon. The spectral parameters are related to these slow-roll 
parameters by [21 [2H] : 

ns = 1 - 6e + 2?7 , 

as = -24e2 + 16er/ - 2^2^ 

13s = -192e^ + 192e^ri - 32eri^ - 24:6^^ + 2ri^^ + 2a^ . (3.3) 



In most (small-fleld) inflation models, eqs.(3.2) imply two strong inequalities between (combi- 
nations of) slow- roll parameters (hierarchy): 

|e| < 1^1, 

\ev\ « \e\- (3.4) 

The flrst relation means that — 1 is essentially determined by r]. Similarly, both relations 
together imply that as is basically flxed by while only the last two terms in the expression 
for (3s cire relevant; these two terms are generically of similar order of magnitude. 

Along with these, another crucial inflationary observable is the influence of gravitational 
waves, relative to density perturbations, on large-angle microwave background anisotropies, 
given by [2] 

-IS-- (-) 
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Combining these equations gives 

as ^ 6(Q' + 4(Q(r^5-l)-2e^ 

-15 Q\ns-1)-2Q {ns-1] 
- {ns - 1)1 + 2a' . 



(3s 



-15 



+ 



as 



(3.6) 
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Figure 2: Contours of q;s + 2^^ (left) and of /3s — 2cr^ (right); the right frame assumes neghgible 
tensor modes, r = 0. 



In the left frame of figure [2] we show contours of 05 + 2^^ in the (ns,r) plane, while contours 
of Ps — 2a' in the {ns,as) plane are shown in the right frame, assuming a negligible tensor- 
to-scalar ratio. We see that the observational constraints on ns and r imply that as + 2^^ is 
very small, roughly —9 x 10~^ < «s + 2^^ < 2.2 x 10~^ if ns and r are within their current 
2 a intervals. Any significant running must therefore be due to [30]. Similarly, —5 x 10"^ < 
Ps — 2a' < —4 X 10~^ if r = and ns and as are within their la intervals. Even using 2a 
intervals and allowing r < 0.17, this range only expands to —0.0053 < Ps ^ 2a' < 0.001, so 
that significant positive running of the running can only be due to a'. 

In slow-roll approximation, the absolute normalization of the power spectrum is given by 



7^c 



(3.7) 



Finally, the number of e-folds of slow-roll inflation that occurred from the time t^, when 
observable CMB scales first crossed the Hubble radius during inflation to the epoch tend when 
inflation ended is given by 

1 /"^* V 

N=^l Y7jd<f>. (3.8) 



"/■end 



v 
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where 0end is defined by max [e((/)end)5 |'7(0end)|] = 1; note that infiation might end through 
dynamics of other fields coupled to the infiaton, as in hybrid infiation. The observationally 
required value of depends logarithmically on the reheating temperature]^ Assuming instan- 
taneous change from infiation to relativistic matter domination, a reasonable range of values of 
the number of e-folds between and tend is taken to be = 54 ± 7 [31]. Requiring baryoge- 
nesis to take place at or above the electroweak scale implies that A^ > 30. A value of A^ ~ 60 
corresponds to a GUT scale reheating]^ 

In the following we study the possibility of PBH formation in two different categories of 
infiation models: small-field models and large-field models. Hybrid models [32] are not studied 
here because in these models, PBH formation can occur by different mechanisms. (For analyses 
of PBH formation in hybrid models, see [33].) As noted above, the spectral index will have to 
increase at very small scales (very large k) in order to allow PBH formation, but we will also 
check whether the models we analyze are compatible with a sizably negative value of as at 
scales probed by the CMB and LSS data, as indicated by eqs.(3.1). 



3.1 Small— field models 

Small-field models are defined as those for which the variation in the infiaton field is less than 
the reduced Planck mass. Typically, e and hence the amplitude of gravitational waves generated 
in such models is small and the spectral index and its running provide the key observational 
discriminators. 

3.1.1 Hilltop/inflection point infiation 

A popular ansatz for the small field infiaton potential is given by [3] 



1 - 



(3.9) 



where Vq, /i and p are positive constantsn This potential is equivalent to the potential V{(f)) = 

A^ — AMp~" — in the literature 134] which can be specialized to several distinct models: e.g. 

n 

hilltop (n = 2 or = 4) and infiection point (n = 3). We consider the case that the dominant 
term is the leading one, Vq. When p is an integer and greater than 2, such a potential may be 
generated by the self-coupling of the infiaton at tree-level. 



For p > 0, the hierarchies (3.4) hold among slow-roll parameters. So the spectral parameters 



^Instantaneous reheating gives the minimum number of e-folds as one looks backwards to the time of 
perturbation production, while a prolonged period of reheating gives a larger number of e-folds. 



^Arbitrarily many e— folds of inflation might have occurred at i < t*, as in "eternal" inflation. N of eq.(3 
is a lower bound on the total number of e— folds of inflation. 

^This potential is unbounded from belove for (j) — > oo. There must be additional terms that prevent this 
Here we follow the usual assumption that these terms do not affect the dynamics of inflation. 
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are given by: 



ns-l 



as 



-2p{p - 1 — 
-2/(p-l)(p-2) 



p-2 



Mp 



-V(p-l)(p-2)2f^^P 



2 (p-2) 



3(p-2) 



/i 



(3.10) 



Inflation ends at 0end ^ and in order to have a small field model we take /x < Mp. Then 



N 



p — 1 



+ 



1 



p — 2 p(p — 2) 



Mp 



p-2 



(3.11) 



For p > 2, the first term in eq.(3.11) can be neglected. We then find, independently of /i: 



ns-l 

as 



p- 


1 2 


p- 


2N 


p- 


1 2 


p- 


2 m 


p- 


1 4 



1 

N 
2 



{ns - 1) 
(ns - 1) 



(3.12) 



p-2m 

In figure |3] the spectral index, its running and its running of running are shown as functions 
of the number of e-folds before the end of inflation, for p = 4. 

It is clear that in this model both ns — 1 and as are negative, but it is not possible to 
reproduce the observed central value of as, which would require as ~ 0{ns — 1) (see figure |3]). 
Moreover, the value of (3s is also negative. So the conclusion is that this model cannot produce 
sufficient high density fluctuations at small scales to produce PBHs. 

In the case at hand, the power spectrum can be calculated exactly as function of A^. Again 



neglecting the first term in eq.(3.11) we have from eq.(3.7): 

1 Vofi^ " 



127r2 J92]Vf6 



'P L 



Np{p 



2)4 



2p-2 
p-2 



(3.13) 



Note that the exponent is positive for p > 2. This implies less power at smaller N, i.e. at 
smaller length scales. 

3.1.2 Running mass inflation 

Another small-field model of interest is the running-mass model [35]. The possibility of PBHs 
formation in this model is studied in detail in [7j. This model is based on the inflationary 
potential 



(3.14) 



where is a real scalar. The potential (3.14) by itself would lead to eternal inflation, so again 



some terms need to be added, and we again assume that these terms do not affect the dynamics 



10 




Figure 3: Illustrating the dependence according to eqs.(3.12) of ns (solid curve), IQas (dashed 
curve) and lOO/S^ (dotted curve), on the number of e-folds before the end of inflation, for the 
fixed value of p = 4. 



of the inflaton during the slow-roll phase. During inflation, the potential is dominated by the 
constant term Vq. Here m^(0) is obtained by integrating a renormalization group equation of 
the form 



din c 



2C ^2 
a m + 

TT 



D 
16^ 



I \ |2 2 
Ay 



(3.15) 



where Pm is the /3— function of the inflaton mass parameter which arises from the gauge inter- 
action with coupling a and from the Yukawa interaction Ay. C and D are positive numbers 
of order one, which depend on the representations of the fields coupling to m is a gaugino 
mass parameter, while is the scalar SUSY breaking mass-squared of the scalar particles 
interacting with the inflaton via Yukawa interaction Ay. 

Over a sufficiently small range of 0, or small inflaton coupling, we can perform a Taylor 
expansion around (p^ 



V = Vo + \ml 



b^ + -c^^ In 



+ '-g<P' In^ 



where 0* is the local extremum of the potential. Here c 



dm"^ 



din I 



(3.16) 

is given by the /3— function. 
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cP"m? 



and g 



is given by the scale dependence of the parameters. By having the 



potential in hand and noting that the hierarchies (3.4) among slow-roll parameters hold in this 
model, we find the spectral parameters, 



ns-l 



as 



L + 1 + f + 3L + 1) 
2c 



2\ 2 



9_ 
2c 



(2L + 3) 



2c 



(3.17) 



cMj 



2\ 3 



2c 



^(3L + 2) + |l(3L^ + 5L + ^ 



where L = In — . Clearly the spectral index is not scale-invariant unless c and g are very close 

(p* 

to zero. 

An important feature of the running-mass model is that because of a consistency relation 
among the spectral index and its running, it is not possible [7] to get large negative running, 
which is presently favored by observation: 

\2 



as > 



{ns 



4 



(3.18) 



On the other hand, we showed in [7] that over a narrow region of allowed parameter space 
this model predicts sufficiently large density perturbations at small scales to allow 10^^ g PBHs 
to saturate the Dark Matter density. 



3.1.3 Inverse power law inflation 

A generic feature of models in nonperturbative gauge dynamics in SUSY [36] is the presence of 



Ar^ 



where the index p and the scale A3 depend on the underlying 



scalar potentials of the form 

gauge group. Like models of hybrid infiation [32], these models are characterized by a potential 
dominated by the constant term Vq and require coupling to another sector to end infiation when 
(f) reaches the critical value (pc Unlike standard hybrid infiation models, models of this type 
postulate a field far from the minimum of the potential. 

We take the potential to be described by a single degree of freedom 0, of the general form 



AP+4 

^^3 



(3.19) 



where the dots represent nonrenormalizable terms suppressed by powers of the Planck mass, 
which are not relevant for the present discussion, but will prevent from "running away" to 
infinity. In the limit -C (0), the term ~ 0~^ dominates the dynamics: 



V(0) 



1 + a 



« (0) 



(3.20) 
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where a 



ATI 



. We assume that the constant Vq dominates the potential, or a ^ (0/Mp)^. 



In this case also hierarchies (3.4) hold among the slow-roll parameters which leads to the 
following spectral parameters: 



71-5 — 1 ~ 2p(p + l)a;^ 



as ^ -2p2(p+ i)(p + 2)a^ 



lis 



/Mf 



2(p+2) 



3(p+2) 



(3.21) 



The number of the e-folds is given by 

N 



p{p + 2) a 



Mf 



p+2 



Mf 



p+2 



(3.22) 



where 0c is the critical value at which inflation ends. The value of 0c is in general determined 
by a coupling of the field to some other sector of the theory which we have here left unspecified. 



Note that from eq.(3.22), for ^ 0c the number of e-folds approaches a constant, which we 
call iVtot, 



tot 



p{p + 2)a \Mp 



P+2 



(3.23) 



This puts an upper limit on the total amount of expansion that takes place during the inflation- 



ary phase, although that upper bound can in principle be very large. Using eqs.(3.21 )-(3.23) 



we can rewrite the cosmological parameters as functions of the number N of e— folds before the 
end of inflation: 



ns-l 



as 



p + 1 



p + 2 



iVtot I 1 



N 



tot 



p + 1 



(3s 



p + 2 
P + l 



tot 



tot 



p + 2{ns-l] 
'p+l 2 



P + 2V (^5-1) 
p+l 



(3.24) 



This model thus predicts ns > I, which is currently disfavored at more than 2 standard 
deviation. Moreover, figure |4] shows that the spectrum becomes scale-invariant towards the end 
of inflation, i.e. it becomes less blue at smaller length scales, as also indicated by the negative 
value of as- In combination with the constraint that \ns — 1| <C 1 at CMB scale, this implies 
that this model cannot accommodate PBH formation. 
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ns 



100 Ps / 



10 as \ 



I 

I 

I 

I 

I 



N 

Figure 4: Spectral parameters as a function of the number of e-folds oc ln(A;) for p = 1. Note 
especially the rapid approach to scale-invariance at short wavelengths (small A^). 



This can also be seen by directly computing the power spectrum as a function on N: 

Vo 



1 2p+2 

'ap{p + 2){N,,,-N)] v+2 



(3.25) 



where we have used (3.22). The power does increase with decreasing A^, but only by a small 



amount. For example, the ratio of the power at the end of inflation (A^ = 0) to that at the 
COBE scale is 

iVtot ; 



COBEJ 



2p+2 
P+2 



(3.26) 



On the other hand, the first eq.(3.24) gives 

p+l 



P + 2 ( Ar,„, 

jVcoBE I -L 



AT, 



COBE 



(3.27) 



Eq.(3.26) can be rewritten as 



tot 



p+2 
g2p+2 _ X _ 



(3.28) 



AT, 



COBE 
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Inserting this into eq.(3.27) finally yields 

p + 2 



p+2 P + 2 

= 1 + 2(p^i) ^coBE [ns{NcoBE) - 1] ~ 2 . 



(3.29) 



The power can therefore only increase by small amount in the course of inflation; in contrast, 
PBH formation would require an increase by a factor 10^ or so. 



3.2 Large— field models 

Large-field models are characterized by the condition | A0| > Mp. Note that a super-Planckian 
field variation is a necessary condition for the generation of an observable tensor-to-scalar ra- 
tio [37]. On the other hand, such large field models raise issues of stability in the presence 
of "quantum gravity" corrections, which are suppressed by inverse powers of Mp. These cor- 
rections should not be important for small-field models, but need not be small for large-field 
models. 



3.2.1 Power— law (a.k.a. chaotic) inflation 

The polynomial potential V(0) = A'' f — ) is equivalent to V((f)) = zz0^ the literature 



. In this model, the hierarchies (3.4) do not hold. We find: 



ns-1 = -p{p + 2) 
as = -2p\p + 2) 



if) 



Ps = -8p'ip + 2) 



3f I ^ 



r = 7p^(^J . (3.30) 

Note that these quantities are independent of the normalization of the potential (described by 
A'^/l^P or, equivalently, by A), but do depend on its shape (described by p) as well as on the 
field value. 

pMp 

Inflation ends at (/)end = — ^ where e = 1. Then, it is straightforward to rewrite the 

v2 

inflation parameters as functions of the number of e— folds, A^: 



2 



77,5-1 



2{p + 2) 
4N + p 



as = -^^^^ = ~^(ns-l? 
' (AN + py p + 2^ ' ^ ' 

_ 64(p + 2) _ 8 

ins-l). (3.31) 



AN + p p + 2 
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Evidently the spectrum is "red" in this model, n^— 1, as and [3s all being negative. However, 
this model also cannot accommodate the current central values, according to which both Ins' — 1| 
and la^l are of order 10^^. Moreover, iVpivot < 50 implies — 1 < 0.040 (0.059) for p = 2 (4), 
i.e. — 1 comes out somewhat below the current central value in this model. 



Computing the power directly from Eq.(3.7), we find: 



1 



127r2 V/iMp 



p + 2 



2p[N 



(3.32) 



This ciecreases quickly towards the end of inflation (A^ 0), again showing that PBH formation 
is not possible in this model. 



3.2.2 Generalized exponential inflation 

We now turn to the generalized exponential potential |39j : 

V{(P) = A^e^"^/^)' , (3.33) 
where p is a positive dimensionless constant and /i is a constant with dimension of mass. In this 



model the hierarchies (3.4) among the slow-roll parameters again do not hold and the values 



of the spectral parameters depend on the field value 



ns-l 



as 



P 



2(p - 1) 



p-2 



P 



p 



3p-4 



{p-2) 



2p-4 



2p\p-l) 



2 



— p(3p — 4) 



4p-6 



2{p - 2f 



3p-6 



2p-2 



(3.34) 



We allow p to be a positive real (not necessarily integer) number. If p is not integer, cf) has to 
be non-negative to get a real potential. In any case the field will roll from larger to smaller 
values during inflation. 



For p > 2, both terms in the first eq.(3.34), or equivalently both e and r/, decrease with 
decreasing cj). The requirements jr^l < 1, e < 1 then yield an upper bound on 0, but inflation 
will never stop once is below this upper bound. This would require an additional mechanism 
to end inflation; we therefore only consider p < 2 here. 

For p < 2, the requirement |?7| < 1 gives a lower bound on 0, which is approximately given 

by 



pIp 



1|M| 



2-p 



(p < 2) . 



(3.35) 



This bound vanishes for p = 1. Eqs.(3.34) show that this choice leads to a constant spectral 



index ns and vanishing as and Ps- This means that inflation does not end for p = 1. Moreover, 
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0.9531 < < 0.9971 , Npi^^, = 45 




0.001 0.002 



0.003 0.004 



0.005 



0.006 



0.007 



Figure 5: Scatter plot of allowed values of as and fis assuming that ns lies in its currently 



allowed 2 a range and 45 e-folds of inflation occurred after the pivot scale for potential (3.33). 



a;^ = is (mildly) in conflict with present data, and a constant in the allowed range will 
not lead to PBH formation. 

If p > 1, the requirement e < 1 implies an upper bound on 0: 



1 



(p>l) 



(3.36) 



The number of e-folds that occur after the inflaton field had a value cj) is given by 



iV(0) 



p{2-p) \Mp 



/i 



2-p 



/i 



2-p 



where 



eq.(3.36) 



is given by (3.35), and for p > 1, (p has to satisfy 



< 0max, with 0max 



(3.37) 



given by 



For p < 1, ns is always less than 1, in accord with observation. The running of the spectral 
index, given by as, is also negative, but we find as > —{ns — 1)^ for all allowed combinations 
of parameters that allow at least 30 e-folds of inflation after the pivot scale. This model can 
therefore not accommodate a sizable and negative value of as, either. 

For 1 < p < 2, ns — 1 can have either sign, while is always positive, in contrast to the 
bound of current data. Quite large and positive as are in principle possible, if parameters are 



chosen such that the two contributions to 77.5 — 1 in eq.(3.34) cancel approximately. However 



such large values of as are definitely in conflict with observation. Moreover, if as > 0.004, 
Ps turns negative, limiting the growth of power at small scales. This is illustrated in figure |5j 
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which shows a scatter plot of allowed values of and assuming that ns lies in its currently 
allowed 2 a range and 45 e-folds of inflation occurred after the pivot scale. 



Using the exact expression (3.7) we numerically find n < 1.15 at scales relevant for the 
formation of 10^^ g PBHs; we saw in section 2 that n > 1.37 is required for the formation of 
such PBHs. This model therefore cannot accommodate PBH formation, either. 

3.2.3 Inflation with negative exponential and Higgs inflation 

Another potential that has been proposed is ^ 

V{(j)) = Vo{l- e-^'^/^*) . (3.38) 
For g > the inflaton field (p rolls towards smaller field values during inflationj^ The 



potential is sufficiently flat only for qcp > Mp, where the hierarchies (3.4) between slow-roll 
parameters hold. The inflationary parameters are 



ns 



-1 ~ -2gV«<^/^^^ 
as -2q'e-"^f/''- = -lins - If ; 

13s ^ -4g6e-39'^/^^p = -{ns- if ; 



where we have approximated the denominators of eqs.(3.2) by Vq; this is appropriate for the 

In Q 

phase of slow-roll where ns — 1, unless ^ 1. Inflation ends at 0cnd = 2Mp . N e— folds 

Q 

before the end of inflation the spectral parameters are given by 

2 

ns — 1 — 



N + 1 
2 



Note that the g-dependence cancels when the spectral parameters are expressed in terms of A^. 

In this model ns — 1 and as are manifestly negative, in agreement with current data. 
However, while ns — 1 also has approximately the right magnitude, \as\ at the pivot scale is 
much smaller than the experimental central value. Moreover, since as and Ps are both negative, 
PBH formation is not possible in this model; this can also be seen from the exact expression 



(3.7), which shows that the power always decreases with decreasing N. 



IOt 



If negative values of <j) are allowed, the potential (3.38) becomes unbounded for (f> ^> — oo. In this case 



additional terms have to be added to the potential, which we again assume to be unimportant during the 
slow-roll phase. 
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A very similar potential describes the Higgs inflation model [lO] where the Higgs boson of 
the Standard Model (SM) plays the role of the infiaton. Starting point of this model is the 
non-minimal coupling of the Higgs field to gravity. The relevant part of the action in the 
Jordan frame is: 



S, 



d Xy 



-9 



2 2 ' 



(3.43) 



where M is some mass parametei , i? is the scalar curvature, h is the Higgs field in the 



unitary gauge and ^ determines the coupling of the Higgs to gravity]^ By making a conformal 
transformation from the Jordan frame to the Einstein frame one can get rid of the non-minimal 
coupling: 



(3.44) 



This transformation induces a non-canonical kinetic energy term for h. It is therefore conve- 
nient to redefine h in terms of the scalar field which casts the kinetic term into the canonical 
form [40] : 



dh 

In terms of this new field, the potential is 

y(0) 



fi2 + Q^-iK^/Ml 



.212 



fi(0)4 4 



(3.45) 



(3.46) 



For small field value, ^ Mp/C,, one has /i ^ and ~ 1; the two frames are indistinguish- 
able so the potential for the field (p is the same as that for the initial Higgs field. However, for 
large values, h > Mp/^/^, one has fi^ ~ ^h'^/M^, and 00] 



h 



Mp 



exp 



(3.47) 



Substituting this into eq.(3.46) we obtain the expression for the potential 



F(0) 



exp 



(3.48) 



Recall that this expression holds only for h ^ Mp / ^/^, which implies that the exponential 
term in eq.( |3.48 ) is small. The square of this term is then even smaller, and can be neglected 
during inflation. The potential therefore effectively almost reduces to the form (3.38 ), wi th 
q = 2/-\/6, except that the exponential term is multiphed with 2 (due to the square in eq.(3.48)). 
This also increases — 1 by a factor of 2: 



ns — 1 — — 



8 Ml 



(3.49) 



^^In the range of ^ of interest to us, M ~ Mp [40] , 

i2jjjggg inflation requires ^ S> 1. This leads to a breakdown of tree-level unitarity at scales well below the 
Planck scale 4T], but according to ref. [H] this does not invalidate the scenario, since the relevant energy scale 
during inflation always remains in the unitary regime. 
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However, this factor of 2 cancels if n5 — 1 is expressed in terms of the number of e— folds of 
inflation that occur after the field had the value 0, i.e. the first eq.(3.40) remains valid in Higgs 
inflation. Moreover, as and (3s are as in eqs.(3.39) when expressed in terms of — 1 or A^. For 
most practical purposes (of inflation), Higgs inflation can therefore be understood as a particle 
physics implementation of negative exponential inflation. 



3.2.4 Natural inflation 

One way to obtain a very flat potential is to consider the natural inflation where a Pseudo 
Nambu-Goldstone Boson (PNGB) is used as inflatonj^ In this model, the inflaton field has a 
particular form of the potential which results from explicit breaking of a shift symmetry: 



1 ± cos a 



f 



(3.50) 



We will take the positive sign in eq.(3.50) and assume that initially (f) <^ f /a. For appropriately 
chosen values of the mass scales, e.g. f/a ^ Mp and A ~ Mqut ~ 10^^ GeV, the PNGB field 
can drive inflation. We set a = 1 for simplicity and treat / as a free parameter. The Slow-roll 
parameters are then given by 



a" 



MpV r sin(0//) y^l/Mp 
T) [i + cos(0//)J ^8VT 
MpV r cos(0//) 1 _ _1 /Mp 
T) Li + cos(0//)J ^"2 VT 
Mp\U sin((/.//) 1^ ^ 1 /Mp 
T) [i + cos(0//)J ^"4VT 
Mpy cos(0//) sin^(0//) ^ 1 /Mp 
/; [l + cos(0//)]^ ^8V/ 



/ 



/ 



(3.51) 



where the approximate equalities hold for (j) <^ f . It is clear that the hierarchies (|3.4|) do not 
hold among the slow-roll parameters 
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and we find the following inflation parameters: 



ns-l 



as 



Mp 

T 

Mp 

T 

Mp 

T 

Mp 

T 



3 - cos(0//) 
1 + cos(0//) ' 
' l-cos(0//) 

[l + cos(0//)]^' 
^ [1 - cos(0//)] [3 - cos(0//)] 

[1 + cos(0//)]=^ 
1 - cos(0//) 



1 + cos(0//) ■ 



(3.52) 



-'^'^Natural inflation can be either a large- or small-field inflation model, depending on the value of /. Here 
we assume / > Mp. 



The flrst strong inequality in (3.4) does hold for </) <C /, but the second one does not even hold in this limit. 
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Inflation ends at \ri\ 
is given by 



1 and the relation between the inflaton field and the number of e-folds 



cos 



/ 



(3.53) 



where y = — e , 

+ 1 



X = 



Mp 
7 



715-1 



as 



/3s 



. Inserting eq.(3.53) into (3.52) yields 



— X 



-Ax' 



-Ax' 



7x' 



2 + y 
2-1/' 

y 



(2-y)2' 

2/(2 + y) 



(2 



(3.54) 



In the course of inflation y increases from a rather small value to y, 



2 + x^ 



end 



at the end of 



1 +X2 

inflation (A^ = 0). Eqs.(3.54) show that ns — l, as and (3s become more negative as y increases, 
indicating that the power is reduced at smaller scales. This can also be seen from the exact 



expression (3.7) which gives 



■Re 



y 



(3.55) 



which decreases with increasing y G [0,2]. PBH formation is therefore not possible in this 
model. 

For fixed N, the spectral parameters are determined by x. \ns — 1| can clearly be made 
as large as desired (with ns < 1) by choosing a large value of x, i.e. a small value of /. On 
the other hand, \as\ reaches a maximum at x^ ~ 2/N for y <^ 1; note that y decreases with 
increasing x^. This gives as ^ —l.b/N"^, i.e. \as\ at the pivot scale cannot be larger than 10~^ 
in this model, well below the current central value. 



3.2.5 Arctan inflation 

Another inflation model which we are interested to study has been introduced in 



1 H — arctan 
vr \ /i 



(3.56) 



This model allows inflation with ~ 1 if /i 3> Mp or ^ /i. However, inflation can be ended 
by the potential (3.56) only if /i < 0.8Mp, since otherwise e, |?7| < 1 V0. A finite period of 
inflation thus requires that ^ initially; at the end of inflation. 



-oo, I.e. 



V ^ 0. 



During the slow-roll phase the hierarchies (3.4) between the slow-roll parameters hold, and we 
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find 



ns-l 



as 



4 /MpV/i 

71 V ' 

12 fMpY ffi 
_72 fMpY fjj 
7 /MpV i^/i 



9 



{ns - If 



(3.57) 



where we have approximated the denominators of eqs.(3.2) by 2Vo, as appropriate for the slow- 
roll phase where 3> /i. In terms of the number N of e-folds of inflation that occurred after 
the inflaton field reached the value 6, we find 



ns 



3N + TI 



(3.58) 



1 /3 

where we have used the fact that infiation ends at 0end — (Mpjj) . This agrees with the 
currently allowed range for 38 < < 95. However, while as is negative, its absolute value is 
only of order 10~^ for allowed values of ns', moreover, since Ps is also negative, PBH formation 



is not possible. Indeed, one can see from the exact expression (3.7) that the power decreases 
steadily during infiation: 



7^c 



127r2M| 



1 + — arctan 



/i 



(3.59) 



This decreases with decreasing for 0end < (f> < oo. 



4 Summary and Conclusions 

In this paper we reviewed the formation of primordial black holes using the Press-Schechter 
formalism. We found that the formation of PBHs with mass larger than 10^^ g, whose lifetime 
exceeds the age of the Universe, will be produced at sufficient abundance to form the cold Dark 
Matter if the spectral index at scale fcpBH is about 1.37 for the threshold value ^th = 1/3. This 
spectral index is much above the value measured at much larger length scales in the CMB. 
PBH formation therefore requires significant positive running of the spectral index when k is 
increased. 

We compared this with the values of the spectral index and its running derived from current 
data on large scale structure. These include analyses of CMB anisotropics from the WMAP 
(7 year) and SPT collaborations, as well as data on baryonic acoustic oscillations and on the 
abundance of clusters, and direct measurements of the Hubble constant Hq. At the pivot scale 
of this data set one finds ns{kpnot) = 0.9751 as central value. The first derivative as{ko) would 
then need to exceed 0.020 if it alone were responsible for the required increase of the spectral 
index; this is more than 3 a above the current central value of this quantity. However, the 
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second derivative (the "running of the running") of the spectral index is currently only very 
weakly constrained. We showed in a model-independent analysis that this easily allows values 
of n(/cpBH) large enough for PBH formation, even if the first derivative of the spectral index is 
negative at CMB scales. 

In section 3 we applied this formalism to a wide class of inflationary models, under the 
constraints imposed by the data mentioned above. We classified the inflation models in small- 
field and large-field models. We have shown that only one small-field model, the running-mass 
model, allows sizable positive running of running of the spectral index, and is thus a good 
candidate for long-lived PBHs formation, albeit only in a narrow range of parameter space. In 
contrast, all the large-field models we studied predict small or negative values for the second 
derivative of the spectral index, and thus predict negligible PBH formation due to the collapse 
of overdense regions seeded during inflation. Recall, however, that other mechanisms of PBH 
formation have been suggested, e.g. in the "waterfall" phase of models of hybrid inflation, or 
during first order phase transitions. 

As a by-product of our analysis, we found that most of the models we studied either predict 
< 1, as indicated by present data, or can at least accommodate it, the single exception being 
inverse power law inflation (a large field model). In contrast, none of the models we analyzed 
allows to reproduce a large negative value of as, as preferred by current data; this confirms 
the results of the general analysis of ref . |15] . If future data confirm with high precision that 
(^s ^ — 0.01, all simple single-field models of inflation would be excluded. Similarly, proving 
conclusively that the second derivative of the spectral index is positive would exclude all the 
large-field models we investigated. Future analyses of the spectrum of primordial density 
perturbations thus hold great promise to discriminate between inflationary scenarios, or even 
to challenge the paradigm of single-field inflation. 
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